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ABSTRACT 

The  advection  of  a  passive  scalar  by  incompressible  turbulence  is  considered  using  recur¬ 
sive  renormalization  group  procedures  in  the  differential  subgrid  shell  thickness  limit.  It  is 
shown  explicitly  that  the  higher  order  nonlinearities  induced  by  the  recursive  renormaliza¬ 
tion  group  procedure  preserve  Galilean  invariance.  Differential  equations,  valid  for  the  entire 
resolvable  wavenumber  k  range,  are  determined  for  the  eddy  viscosity  and  eddy  diffusivity 
coefficients  and  it  is  shown  that  higher  order  nonlinearities  do  not  contribute  as  k  — »  0,  but 
have  an  essential  role  as  k  kc,  the  cutoff  wavenumber  separating  the  resolvable  scales  from 
the  subgrid  scales.  The  recursive  renormalization  transport  coefficients  and  the  associated 
eddy  Prandtl  number  are  in  good  agreement  with  the  A:-dependent  transport  coefficients 
derived  from  closure  theories  and  experiments. 
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1  Introduction 


The  turbulent,  transport  of  a  passive  scalar,  while  serving  .sourni  /K'dagogical  purposes,  is 
also  of  interest  in  the  spreading  of  temperature,  humidity  and  jjollution  in  the  atmos]>here 
as  well  as  iji  other  problems  {Csanady' ).  Here  we  shall  apply  recursive  renormalization  group 
(RNC)  procedures  to  the  subgrid  modeling  of  a  pa.ssive  scalar  field  7'(k. /)  being  advected 
by  a  turbulent  Navier-Stokes  velocity  field  u(k,  i).  Subgrid  modeling  is  necessary  for  tlu‘ 
high-Reynolds  number  turbulent  flows  of  interest  becau.se  of  the  limitations  of  current  and 
foreseeable  supercomputers^.  Another  advantage  of  considering  the  problem  of  passive  scalar 
transport  is  that  the  spectral  transport  coefficients  (eddy  diffusivity  and  eddy  visco.sity) 
de+ermined  f’-om  our  RNG  theory  can  be  compareci  to  those  arising  from  closure- based 
theories^’”*.  It  should  be  noted  that  the  trans])ort  coefficients  in  these  closure  theories  are 
determined  over  the  whole  resolvable  scales. 

Recently,  two  distinct  approaches  of  RNG  to  fluid  turbulence  have  ari.sen:  one  leased  on 
the  work  of  Forster  et.  al.^,  called  c-RNG,  and  the  other  based  on  Rose*’,  called  recursive'- 
RNG.  Some  aspects  of  these  two  approaches  have  been  discussed^.  In  particular,  we  point 
out  here  that  in  e-RNG,  a  small  parameter  e  is  introduced  through  the  forcing  correlation 
function.  Yakhot  &  Orszag®  have  extrapolated  c  1  to  e  — >  4  in  order  to  reproduce 
the  Kolmogorov  energy  spectrum.  Furthermore,  it  is  also  necessary  to  take  the  distant 
interaction  lirndt®,  0.  Thus,  it  is  difficult  to  compare  the  transport  coefficients  generated 
by  Kraichnan®  and  Chollet'*,  with  that  determined  from  c-RNG. 

In  this  paper,  we  continue  our  application  of  recursive  RNG^*^"”  to  turbulence.  The 
basic  differences  between  the  two  RNG  procedures  are  that  in  recursive  RNG: 

(i)  The  e-expansion  is  not  applied. 

(ii)  The  turbulent  transport  coefficients  are  determined  for  the  whole  resolvable  wavenum¬ 
ber  scales, 

(iii)  Higher  order  nonlinearities  are  generated  in  the  renormalized  momentum  equation. 

In  Sec.  2  we  derive  the  renormalized  evolution  equations  for  the  passive  scalar  7’(k, /) 
and  the  fluid  velocity  u(k,  t)  as  well  as  the  recursion  relations  from  which  the  eddy  diffusivity 
and  eddy  viscosity  can  be  determined.  Because  of  the  presence  of  higher  order  nonlinearities 
in  the  renormalized  equations,  it  is  not  apparent  that  Galilean  invariance  is  still  preserved. 
These  questions  are  addressed  in  Sec.  3  where  we  prove  that  the  RNG-evolution  equations  are 
Galilean  invariant:  a  property  deemed  necessary  in  any  subgrid  model'^.  It  has  bc’en  found 
very  difficult  to  find  fixed  points  for  the  RNG-difference  recursion  relations  if  the  subgrid 


shell  thickness  is  chosen  too  small'®’*.  If  recursion  RNC  procedures  are  to  he  em|)loyed 
successfully  in  more  complicated  flow  problems,  then  it  is  necessary  that  these  dstference 
recursion  relations  be  simplified.  In  Sec.  4  we  proceed  to  the  differeistial  limit  of  these 
recursion  relations,  paying  careful  attention  to  the  k  Q  limit.  We  show  that  the  iiigluu 
order  RNG-induced  nonlinearities  do  7iot  contribute  to  the  k  0  limit  of  the  RNCl  rc'cur.sion 
relations,  but  play  a  significant  role  for  k  — *•  kc,  where  kc  is  the  wavenumber  that  separates  the 
resolvable  scale  from  the  subgrid  scale.  In  the  Appendix,  we  contrast  our  ordinary  differential 
equations  for  the  RNG  eddy  viscosity  and  diffusivity  with  that  generated  by  Yakhot-Orszag'^ 
by  their  c-RNG,  a  theory  that  is  valid  only  in  the  distant  interaction  limit  k  — >  0.  In  contrast, 
in  the  recursive  RNG  approach,  the  differential  equations  for  the  eddy  transport  coefficients 
are  valid  over  the  whole  resolvable  scales  and  not  just  at  k  0,  and  no  (  expansion  is  needed. 
The  turbulent  transport  coefficients  for  the  second  moments  (i.e.,  for  the  time  evolution  of 
6a0(k,t)  — <  itQ{k,  k,  i)  >  and  the  scalar  variance  0(k,t)  =<  r(k,  k,  t)  >  ) 

are  determined  in  Sec.  5.  In  particular,  the  importance  of  the  RNG-induced  higher  order 
nonlinearities  is  very  apparent.  The  spectral  eddy  viscosity,  diffusivity  and  Prandt!  number 
are  derived  in  Sec.  6,  while  in  Sec.  7  we  present  our  conclusions. 

2  Renormalized  momentum  equation  for  velocity  and 
passive  scalar 

We  consider  a  passive  scalar  ^(k,^)  being  advected  by  incompressible  turbulence 


[-^  + fiok‘^]T{k,t)  = -ike  J (1) 

with  the  turbulent  velocity  field  u(k,  f)  being  determined  from  the  Navier-Stokes  equation 

[^  +  i^ok'^]ua{k,t)  -  Ma0-y{k)  j  d^jup{},t)u^{k- 3,  t) -jr  fe{k,t).  (2) 

Summation  over  repeated  subscripts  is  understood,  and 

^oc0'f{k^  kpDe'yi^k^  T  k^De^^k'^^  and  D^p^^k^  —  k^k^jk  .  ^3) 

Here  /^o  is  the  molecular  diffusivity,  vq  the  molecular  viscosity  and  /<,  a  random  forcing  term. 
The  forcing  correlation  is  given  by 

<  fe{k,  t)f0{k\  t')  >=  Dok-y  Dep{k)8{k  +  k')S[t  -  t')  (4) 


where  Do  denotes  the  intensity  of  the  forcing,  and  y  is  an  appropriately  chosen  exponent  so 
as  to  recover  the  Kolmogorov  energy  scaling  in  the  inertial  range  (y  =  3). 
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2.1  An  outline  of  the  recursive  RNG  procedure 

Since  the  details  of  the  recursive  RNG  procedure  for  Navier-Stokes  turbulence  have  been 
presented  before®’*®'",  we  just  briefly  outline  the  steps  here; 

(i)  The  subgrid  wavenumber  region  {kc,kd)  is  partitioned  into  A'^-shells 

kc  =  I’/v  <  fciv-i  <  ••••  <  kx  <  ko  =  kd  (5) 

where  kc  is  the  wavenumber  separating  the  resolvable  from  the  subgrid  scales  and  kd 
is  the  Kolmogorov  dissipation  wavenumber.  A:„  =  n  =  0, N ,  where  /  is  a 

factor,  0  <  /  <  1,  measuriu;:^  the  coarseness  of  the  subgrid  partitioning.  The  limit 
/  ^  1  corresponding  to  a  differential  partitioning  of  the  subgrid  region  {N  — ♦  oo). 

(ii)  The  subgrid  modes  for  the  first  shell,  ki  <  k  <  k^,  are  eliminated  from  the  resolvable 
scale  equation  by  the  solution  of  the  subgrid  scale  equation. 

(iii)  A  subgrid  scale  average  is  performed  over  the  resultant  resolvable  scale  equation.  This 
will  result  not  only  in  the  introduction  of  the  subgrid  scale  energy  (or  equivalently, 
forcing)  spectrum,  but  it  also  results  in  a  new  triple  nonlinearity  and  nonlocal  eddy 
damping  function  in  the  resolvable  momentum  equation  {k  <  kc)- 

(iv)  The  above  steps  are  repeated  for  each  successive  subgrid  shell  until  all  the  subgrid 
scales  have  been  removed. 

(v)  Since  the  subgrid  scales  evolve  on  a  faster  time  scale  than  the  resolvable  scales,  a 
multiple  time-scale  analysis  can  be  performed  to  simplify  the  eddy  damping  function. 
The  resultant  eddy  viscosity  is  a  fixed  point  of  an  integro-difference  recursion  relation. 

(vi)  The  recursion  relation  for  the  eddy  viscosity  and  the  renormalized  Navier-  Stokes  equa¬ 
tion  are  rescaled. 

It  should  be  emphasized  that  there  are  two  singular  limits^:  /  — +  1  and  ^  +  0.  A  careful 
analysis  must  be  done  regarding  these  two  limits  and  the  associated  averaging  operations. 
We  will  address  this  issue  here  in  the  present  paper. 

2.2  An  asymmetry  in  the  renormalized  passive  scalar  equation 

The  details  of  the  implementation  of  the  recursive  RNG  procedure  to  the  advection  of  a 
passive  scalar  is  a  straightforward  generalization  of  that  for  Navier-Stokes  turbulence  (see 
e.g.  Hossain*®)  and  so  will  not  be  presented  here.  Here  we  comment  on  a  symmetrization 


procedure  that  is  standard  when  dealing  with  Navier-Stokes  turbulence,  but  which  can  not 
be  applied  to  the  passive  scalar  problem. 

Consider  the  removal  of  the  first  subgrid  shell  and  use  the  usual  notation 


_  r  u^(k,/)  if  A-,  <  A  <  Ao  ; 
l«,f(k,  0  ifA<Ai. 


Tile  f  \  —  i  ^^(k-.0  if  ^‘i  <  A  <  Ao  ; 

^  *  ir<(k,t)  ifA<A,.  ' 

We  find  for  A  <  A’j,  the  resolvable  scale  passive  scalar  and  Navier-Stokes  equation  can  be 
written  as 

+«;(k-j,0r>(j.0  +  Kj(k-j,()7’'’(j.()l.  (8) 

and 

^  +  ^'oA^]u<(k,<)  =  f^{kj)  +  Mc0^{k)  J j[up{}J)u'^{k  -  iU) 

-f2u|(j,t)w<(k  -  u^(j,0w^(k  -  j,0]-  (9) 

The  factor  2  in  the  Navier-Stokes  Eq.  (9)  arises  from  the  symmetry  in  the  j  k  —  j 
interchange. 

We  assume  isotropy  for  both  the  velocity  field  and  passive  scalar,  so  that  the  subgrid 
velocity  -  passive  scalar  correlations  are  zero  (Lesieur*"*): 


<  «>7’>  >=  0, 


where  <  . . .  >  represents  averaging  over  the  subgrid  scales. 

On  applying  the  recursive  RNG  procedure  to  eliminate  the  subgrid  fields  in  Eqs.  (8) 
and  (9),  the  second  term  on  the  RHS  of  Eq.  (8)  requires  special  attention.  In  Ref.  6,  the 
frozen  flow  velocity  field  is  prescribed,  so  that  this  term  plays  no  role  in  the  renormalization 
procedure.  However,  here  we  are  considering  a  passive  scalar  field  being  advected  by  a 
turbulent  velocity  field  which  is  itself  determined  from  Navier-Stokes  turbulence.  As  a  result, 
this  term  must  be  treated  on  the  same  footing  as  the  others.  For  this  second  term  on  the 
RHS  of  Eq.  (8)  we  have  the  wavenumber  restrictions,  j,  A  <  Ai  while  k\  <  lk-j|  < 
Aq  is  in  the  subgrid  scale.  This  can  only  be  achieved  when  A  and  j  are  located  near  Aq. 
Thus  u^(k  —  j, t)T'''{j, <)  with  the  wavenumber  restriction  discussed  above  is  different  from 
0^  ^(k  —  j,  t)  with  A,  jk  —  jl  <  A]  and  A|  <  j  <  Ao.  Hence  one  can  not  interchange  the 
variables  j  and  k  -  j  by  standard  .symmetry  arguments  as  done  by  Hossaiid^. 
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2.3  The  renormalized  Navier-Stokes  and  passive  scalar  equations 

After  removing  the  subgrid  shell,  the  renormalized  passive  scalar  equation  takes  the  form 


+  ^„(i)e]7-(k,0  =  -a,,/ d^jn<{k  -  j.i)7'<(k,i) 

/ d=j,ej'  -j',or<(j',() 

*'  Mn~/i  w  )J 

where  the  last  term  on  the  RHS  differs  from  that  of  Hossain’^. 


The  restriction  on  the  wavenumbers  are  the  following: 

<  j  <  ^>-1  in  the  second  integral 
k;^  <  [k  —  j|  <  in  the  third  integral. 

The  other  wavenumber  constraints  are  as  indicated  by  the  superscript  on  the  fields  u  and  7 . 
The  renormalized  Navier-Stokes  equation  has  the  following  form”; 

[a/at  +  t',(J.-)t*]u„(k,  ()  =  /<,(k,i)  +  M„fl,(k)  J  <pju^(j,t)u^{k-i,t) 
+2M,«,(k)^  / - ^MMvO)“J(j',<)<.(j-j'.<X(k-j.')  (12) 

/i=l  ^n-h\J)3 

where  j  is  restricted  to  the  subgrid  shell  in  the  second  integral.  Again,  all  other  wavenumber 
constraints  are  as  indicated  by  the  superscript. 

2.4  Recursion  relations  for  eddy  viscosity  and  diffusivity 

Although  the  second  term  on  the  RHS  of  Eq.  (8)  contributes  a  new  triple  nonlinear  term, 
it  does  not  make  any  contribution  to  the  renormalized  eddy  diffusivity  in  the  monrentum 
equation  in  the  process  of  removing  the  next  subgrid  shell.  The  reason  is  the  following: 

2"'^  term  ~  j)  <  0<(k  -  j  -  j',  <)  >  T<{lt)  ^  0  (13) 

since  the  ensemble  average  will  generate  a  delta  function  ^(k  — j)  while  k— j  is  in  the  subgrid 
range.  This  is  impossible,  and  so  this  second  term  can  not  contribute  to  the  eddy  diffusivity. 

After  the  removal  of  the  (n  -f  l)‘^  subgrid  shell,  the  spectral  eddy  viscosity  in  the  renor¬ 
malized  momentum  equation  is  determined  by  the  recursion  relation” 


where 


and 


=  i/M  +  Sun{k) 


—  f  ffi  ■  jl _ 


(14) 


(15) 


L{k,j,q) 


kjjl-  2^)[2:q^  -  kj] 


(16) 


with  k  •  j  =  kjz.  This  difference  equation,  after  rescaling,  has  been  solved  by  Zhou  et  al." 
and  fixed  points  were  readily  determined  for  finite  /  <0.7.  However,  it  was  very  difficult  to 
determine  fixed  points  for  finer  subgrid  partition  factor  /  >0.7.  In  Sec.  4  we  shall  pass  to 
the  differential  subgrid  limit  /  — »  1  and  determine  an  ordinary  differential  equation  (o.d.e) 
for  the  renormalized  eddy  viscosity  over  the  entire  resolvable  scale  which  can  be  readily 
integrated. 

In  a  similar  fashion,  the  spectral  eddy  diffusivity  in  the  renormalized  passive  scalar  equa¬ 
tion  can  be  shown  to  be  given  after  the  removal  of  the  (n  -f  1)‘^  subgrid  shell,  by 


fj'n+lik)  =  fJtnik)  +  Sfln{k)  (17) 

where 

The  renormalized  eddy  viscosity  and  diffusivity  are  defined  as  the  fixed  point  of  these 
recursion  relations. 

2.5  Rescaling  of  the  recursion  relation  and  momentum  equations 

From  the  self-similarity  properties  of  the  forcing  and  energy  spectrum  in  the  subgrid  range, 
we  expect  that  the  viscosity  i/^+x  to  be  simply  related  to  i/„  for  large  n,  while  the  diffusivity 
is  simply  related  to  A  rescaling  can  be  performed  on  the  recursion  relation.  In 
particular,  consider 


and  define 


k  *  ^n+l k 

f„(i)  =  for  i  <\ 


(19) 
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Unless  mentioned  otherwise,  we  drop  the  tilde  notation  on  the  wavenumber,  and  note 
that  0  <  /:  <  1.  The  recursion  relation  for  the  eddy  viscosity  becomes 


where 


Shim  =  / <^j - , 

£l  [  •  r-My+l)/3 _ ~  jl  ^ _  (22) 

and  the  summation  term  arises  from  the  triple  nonlinearity  induced  by  the  recursive  RNG 
procedure.  The  recursion  relation  for  the  eddy  diffusivity  is 

Mk)  =  +  «A.(/t)l  (23) 


s,Mk)  =  - -477^3 - 

k  i  ‘  ' 

where  again  the  second  term  on  RHS  of  Eq.  (24)  arises  from  the  induced  triple  nonlinearities. 
We  have  introduced  the  parameter  m  =  5/3.  These  equations  are  valid  for  any  k  in  the 
resolvable  scales:  Q  <  k  <  kr.  In  the  limit  k  —*  0,  the  triple  term  contribution  — >  0,  as  will 
be  shown  in  Sec.  4.1. 

The  final  renormalized  passive  scalar  equation  is 


j^+f,ik)k^\T(k,t)  =  -mj  <P]u<(k-i,t)T<(k,t) 


li(K)h'  ■'  J  > 

. -sa.  I  -  j  -  j',Or"(j.O 


while  the  final  renormalized  Navier-Stokes  equation  is 


i/{kc)kc 


[d/dt  +  u{k)k'^]ua{K  t)  =  /a(k,  t)  +  M„0^{k)  J  (Pju^  (j,  t)u^{k  -  j,  t) 

■i-75rMa0-7(k)  I d^j<^j'-:^Mf}0>Y{j)u^,{y,t)u<,{}~2\t)u<{k-j,t). 
.)kc^  J  ^ 
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3  Galilean  invariance  of  the  renormalized  Navier-Stokes 
and  passive  scalar  equations 

wo  proceed  further,  we  turn  tuir  attention  to  the  (nu'stioii  of  th<‘  (lalilean  iii\ai  iaii<  <■ 
of  thf'  renormalized  Navier-Stoke.s  ami  passive  scalar  ('(juations  (2o)  and  {'Hi).  The  inijn)! 
tance  of  (lal'U'an  invariance  in  turbulence  modelling  has  been  en>|)hasized  by  Speziale''^, 
do  be  consistent  with  the  l>asic  physics,  it  is  recpiired  that  the  descri|)fion  of  the  turlm- 
lence  be  the  same  in  all  inertial  frames  of  referenc<‘.  The  appearanct'  of  the  trij)le  mjiilijiear 
term,  which  is  a  function  of  the  rt'solvable  scales  velocity  fields,  makes  the  property  of  the 
(ialilean  invariance  of  our  recursive  HN(I  procedure  not  transparent.  VVe  now  show  that 
both  the  renormalized  Navier  -Stokes  e(|uation  and  the  renormalized  |)assive  scalar  e(|uation 
are  (ialilean  invariant. 

3.1  Galilean  invariance  in  Navier-Stokes  equation;  A  review 

The  CJalilean  transformai icui  is 


/  ->  r 


Thus,  one  has 


()  d  '() 


.  _  -  '<)  d  ()(),,  <) 

u-u  -u„. 

While  the  (Jalilean  transformation  for  the  Navier-.Stokes  e(|uation  in  physical  space  is 
trivial,  the  (Jalilean  transformation  in  wavenumber  space  is  less  obvious,  due  to  the  lack  of 
differential  operations.  For  convenience,  we  first  review  how  (Jalilean  invariance  is  preserved 
for  the  Navier-Stokes  equation  in  the  W'avenumber  spac*-. 

Under  th<’  (Jalilean  transformation,  the  LHS  of  the  Navier-Stokes  e<|uation  [cf.  Fcp  (2)] 
becomes 


df 


=  -f-  /y,r^q,(k*,0 


vvliert'  in  the  last  step,  we  hav('  used  the  th<'  b  function  |)roperty  A’^<S(k*)  =  0. 

Under  the  (lalilean  transformation,  the  RHS  t)f  the  Navier-Stokes  equation  [c.f.  Eq.  2] 
Ix'comes 


8 


J  ~  J'-o  -  -  j*)] 


(30) 


wliere  vve  have  used  the  property  of  the  6  function,  the  incompressible  condition,  and 


A/a/i^(^-‘)fw';(k*,0  =  Uo,3k:;uJk\t)l{2t) 


(31 


Thus,  as  expected,  the  Navier-,Stokes  equation  is  invariant  under  a  (laliU’an  transforma' 
tion  due  to  the  cancellation  of  the  second  term  on  the  RHS  of  Eqs.  (29)-(30). 

3.2  The  renormaii7ed[  Navier-Stokes  equation  under  a  Galilean 
transformat  ion 

To  show  that  the  renormalized  Navier-.Stokes  equation  is  inv'ariant  under  a  Galilean  trans¬ 
formation,  we  need  only  consider  the  recursive  RNG  induced  triple  nonlinear  term,  denoted 
by  NSr. 

A'6Vs2M.«,(fc)  J  (32) 

It  is  important  to  note  that  j  is  in  the  subgrid. 

Under  a  Galilean  transformation,  Eq.  (32)  becomes 

■  i*/j  )| 


NS-r  -  2M,0,{kn  I  -  j'*,0  -  ^W'<*^(j*  -  j")] 


K(r,o  -  cw<^(jnK(k* -y,t)  -  uoySik^-n 


(33) 


Since  j"  is  in  the  subgrid  scale,  while  /*  and  k'  are  in  the  supergrid,  6(k*  —  j*)  and 
(5(j*  ~  j'*)  can  never  be  simultaneously  satisfied.  As  a  result. 


Nsy  =  2M„,„(e)  /  <f'r<pnff^^^f^uAy  -r.o 

K.(r,o-Co,'<(r)i<(k--j-,n 


(:M) 
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Now  only  one  term  in  Eq.  (34)  could  violate  the  Galilean  invariance  of  the  renormalized 
Navier- Stokes  equation.  However, 


j - r.t)  ^ 


135) 


This  is  not  permissible  since  Ug,  =  u^,’  and  j*  is  restricted  to  the  subgrid,  riius  NSr 
NS^.  Hence  the  triple  term  is  Galilean  invariant. 


3.3  Galilean  invariance  in  the  renormalized  passive  scalar  equa¬ 
tion 

The  renormalized  passive  .scalar  equation  has  two  triple  nonlinear  terms.  The  proof  of  tl)e 
Galilean  invariance  of  the  renormalized  passive  scalar  equation  ]}roc<H'ds  in  a  similar  manner 
to  that  for  the  renormalized  Navier-Stokes  equation. 

For  the  first  triple  nonlinear  term,  labelled  PSt\,  after  a  Galilean  transformation,  we 
found 


p$r,  ~  j  d^rd^riu<-{k-  -r.()  +  OMiu-  - 

=  /rf=;V=ur(k-  - 


since  is  in  the  subgrid  while  are  in  the  resolvable  scale.  Thus  the  6(k’  —  j*)  and 

^(j*  —  )  can  never  be  satisfied  simultaneously. 

The  second  triple  nonlinear  tcm  has  the  following  structure  after  the  Galilean  transfor¬ 
mation 


psf,  ~  j  ^r^riura'-j) + co<,i(r)iK-(k-  -  j- -y.t) + ucav  -  j-  -  j'-)ir<(j.o 

^  J  rf=j-d»[uy{j'-)  +  i/„;,^(j'-)i[uj-(k-  -  j-  - 

-j-  -r,f)r<o.<)(37) 


where  the  last  two  steps  follow  from  the  wavenumber  constraints,  are  in  the  resolv¬ 

able  scales  while  |k— j|  is  in  the  subgrid  scale.  Specifically,  the  first  step  follows  since 
(J(k*— j*  —  j'*)  can  never  be  satisfied.  The  second  step  follows  since  6(j'*)  would  force 
u^*(k*  —  j*  —  j'*)  — +  u;^*(k*— j*).  This  is  not  permissible  since  jk  — jj  is  in  the  subgrid 
while  u*,  by  definition,  is  in  the  resolvable  scale.  Thus,  the  renormalized  passive  scalar 
equation  is  also  Galilean  invariant. 
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4  Differential  equations  for  the  renormalized  eddy 
viscosity  and  diffusivity 

The  differential  limit,  /  — »  1,  is  singular  and  hcis  been  discussed  recently^.  In  particular,  it  is 
related  to  the  assumption  of  local  versus  non-local  interactions  in  k.  In  this  section  we  will 
calculate  the  eddy  viscosity  and  diffusivity  under  the  differential  ecpiation  limit  for  recursive 

RNG. 

For  recursive  RNG  we  will  find  that  the  differential  equations  hold  throughout  the  resolv¬ 
able  wavenumber  range  0  <  This  should  be  contrasted  with  f  -  RNG  eddy  viscosity 

differential  equation  which  is  valid  only  in  the  k  0  limit®. 

4.1  The  distance  interaction  approximation,  k  0 

Consider  the  resolvable  scale  Navier-Stokes  equation,  Eq.  (9), 


u<(k,t)  =  Mapy{k)  j  -  j,0 


The  first  and  third  terms  on  the  RHS  of  (9)  are  symmetric  in  j  and  |k  —  j]  in  terms  of 
their  respective  wavenumber  constraints  in  wavenumbers.  As  a  result,  the  di.stant  interaction 
limit  fc  — >  0  has  no  effect  on  the  existence  of  these  terms  which  will  give  rise  to  the  standard 
quadratic  nonlinearity  and  eddy  viscosity,  respectively.  However,  the  second  term  on  the 
RHS  of  (9)  has  the  following  constraint:  j  is  in  the  subgrid  while  |k  —  j|  is  in  the  resolvable 
scales.  Specifically,  the  consistency  condition  requires  that,  for  small  j  satisfies 


j  >  kc  and  j  <  K.  +  kz.  (38) 

Since  lz\  <  1,  the  range  of  integration  must  be  0{k). 

Thus,  the  second  term  on  the  RHS  of  Eq.  (9)  can  not  contribute  in  the  limit  — »  0  since 
the  integrand  is  bounded.  A  similar  conclusion  can  be  drawn  for  the  second  and  third  terms 
in  the  renormalized  passive  scalar  equation.  Hence,  the  triple  terms  will  not  contribute  to 
the  renormalized  momentum  equations  and  recursion  relation  for  the  transport  coefficients 
in  the  distant  interaction  limit.  A:  — >  0.  However,  they  will  contribute  to  the  renormalized 
Navier-Stokes  and  passive  scalar  equations  for  0  <  A:  <  Ar^. 
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4.2  Test  of  the  conclusion  in  §4.1  from  numerical  simulation 
databases 

The  conclusions  of  the  last  subsection  can  be  tested  directly  using  numerical  simulation 
databases.  Indeed,  energy  transfer  and  eddy  viscosity  can  be  analysed  using  results  from 
numerical  simulations  by  introducing  an  artificial  cut  at  a  wavenumber  k\.  that  is  smalh'r 
than  the  maximum  resolved  wavenumber  k,„  of  the  simulation.  With  this  fictitious  separation 
between  the  subgrid  and  re.solvable  scales,  it  is  possible  to  evaluate  the  effect  of  the  su!)grid 
kc  <  k  <  kjn  on  the  resolvable  scales  To  facilitate  comparison  with  tlie  recursive  UNCI 
analysis,  we  consider  separately  the  contribution  to  the  energy  transfer  and  eddy  viscosity 
from  the  second  and  third  term  on  the  RHS  of  Eq.  9.  We  form  an  energy  e<piation  from 
the  momentum  equation  and  introduce  tlie  following  notation:  {k)  is  the  spectrum  of 

energy  transfer  to  mode  k  resulting  from  interactions  between  modes  with  wavenumbers  less 
than  kc\  T^^(k)  and  T^^{k)  represent  similar  contributions  from  interactions  with  otir  or 
both  modes  above  the  cutoff  k^  respectively.  The  equivalent  contributions  to  eddy  vi.scosity 
in  the  energy  equations  are  v^^{k)  =  -~T^^{k)f2k^ E(k),  and  i/^^{k)  =  {k)/2k^ E{k). 

To  determine  the  behavior  of  the  energy  transfer  and  eddy  vi.scosity  i/^^(k)  and  i/^^{k) 
we  measured  them  in  flow  fields  obtained  from  numerical  simulations  on  128^  meshes  of 
forced  turbulence.  The  forced  flow  dataset  was  generated  by  C'hasnov’'*  in  a  large-eddy 
simulation  (LES)  of  the  Kolmogorov  inertial  range,  using  a  subgrid  model  derived  from  the 
stochastic  equation  that  is  consistent  with  Eddy-damped-quasinormal  Markovian  (EDQNM) 
approximation. 

In  Fig.l,  we  present  a  numerical  measurement  of  o^'^{k).  It  demonstrates  that  the  second 
term  on  the  RHS  of  Eq.  (9)  does  not  contribute  to  the  energy  transfer  process  as  k  — »  0, 
consistent  with  our  analysis. 

4.3  The  differential  equation  limit,  /  — ^  1 

We  now  derive  the  differential  equation  for  the  transport  coefficients  for  finite  k,  0  <  k  <  U 
where  the  wavenumbers  are  normalized  with  respect  to  the  cutoff  wavenumber  k^  =  I.  The 
o.d.e.  in  the  distant  interaction  limit  (k  =  0)  will  be  derived  in  the  next  subsection.  After 
the  rescaling,  we  rewrite  the  recursion  relation  in  the  form 

(39) 

For  /  — >  1,  the  number  of  iterations  n  — +  oo.  Similarity  considerations  lead  to 

— +  i'{k),  n  — >  oo.  (40) 
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Let  A  =  I  —  /.  the  LHS  of  Eq.  (39)  becomes 


//(A-)  -  [1  -  -  A)] 

As  noted  earlier^,  the  partial  average  of  Rose**  must  be  employed  in  order  to  insure  the 
existence  of  the  differential  limit.  The  partial  average  is  introduced  since  the  distinction 
between  the  resolvable  and  subgrid  scales  become  fuzzy  in  the  limit  of  a  differential  subgrid 
partitioning,  f  —*  1. 

Following  Rose'’,  we  first  change  the  variable  from  j,z  to  j,q  =  jk— jj  ,  with  djdz  ~ 
{qlkj)djdq,  so  that  the  RHS  of  Eq.  (39)  becomes 


where 


6„{k)  _*j|)]k  -  j|2|k  -j|» 

I  /  dida(  —  ) _ _ /^\(j,+i)/3 

Ji<q<i+k  k^iy'^{\)q^~^  Ji<j<i+k  ^  ^ 


L{k,\,q)  =  k(\-z^){k-zq^)/q^ 


L{k,j,\)^kj{\-z^){kj-z).  (44) 

Here  one  has  set  the  coefficient  Dq  =  'IttDo  =  1  (Zhou  et  al.”). 

As  a  result,  the  fixed  point  renormalized  eddy  viscosity  i/{k)  is  determined  from  the  o.d.e. 
at  0(A) 

+  B,(k)]  (45) 


where 


BAk)=^j^*‘‘djHk.jJ)fi- 


(y-2;/3 


Here,  2  is  evaluated  at  j  =  1  and  <7=1,  respectively  in  the  L{k,j,q)  expression. 
The  fixed  point  o.d.e.  for  the  eddy  diffusivity  is  that  given  by 


(48) 


where 


+  t±la{k)  =  — 
dfc  ^  3  /x(l) 


K(/t)  +  5,(A:)j 


A, 


<<i<i+*,- 


<i<7 


sin^(^,  q,  1) 


,m+l 


(49) 


«*•(*)  =kf  ‘‘i  (50) 

and  sin^(^,j,  g)  is  the  square  of  the  sine  of  the  angle  defined  by  the  k  and  q  legs  of  the 
k,j,q  wave-vector  triangle.  Note  that  Eqs.  (48)-(50)  are  identical  with  the  o.d.e.  of  eddy 
diffusivity  which  Rose®  derived  for  a  prescribed  frozen  velocity  field. 


4.4  Differential  equations  in  the  fc  — v  0  limit 

In  the  A;  — »  0,  we  have  seen  that  the  triple  term  does  not  contribute  to  the  eddy  viscosity.  As 
a  result,  the  recursion  relation  will  now  contain  only  the  usual  quadratic  contribution.  We 
further  simply  the  analysis  by  taking  the  standard  subgrid  linear  propagator"  —  jj)  = 

The  limits  of  the  integration  are  given  by 


l<j~-kz<\ff,  \  ■+ kz  <  j  <  Iff -k  kz  (51) 

Thus,  the  RHS  of  Eq.  (39)  becomes 


Sunik)  =  M-A-B  (52) 


where  the  integral  limits  for  these  terms  are 

/■!//  /■! 

1 

for 

M 

(53) 

fl+kz 

/„  '*■’ 

for 

A 

(54) 

r  dz  dj 

for 

B. 

(55) 

J-l  Jilf+kz 

Terms  A  and  B  are  the  corrections  to  the  symmetric  term  M.  They  are  important  for 
a  finite  bandwidth  /.  However,  it  is  easy  to  show  that  A  +  5  =  0  for  /  — »  1  in  the  fc  0 
limit.  Hence 


/_(  -  1)]  =  A 
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15i/2(l) 


(56) 


14 


while  the  LHS  of  Eq.  (39)  yields 


dk 

since  is  bounded  as  k 
Thus,  as  A:  — >  0, 


jdp{k)  ,  3/  +  1  ,,,  y  +  1 
k  r:  I - -  i>{k)  -  t'ik), 


as 


0. 


u{k  — >  0)  = 


3  8 


(57) 


y  +  1  15i/2(l)- 
Again,  one  sets  the  coefficient  =  2t:Dq  =  1. 

A  similar  anlysis  can  be  performed  on  the  fixed  point  o.d.e  for  the  eddy  diffusivity.  This 
was  not  performed  by  Rose®  who  did  not  consider  the  A:  — ♦  0  limit  carefully.  Again,  as  A-  -+  0, 
the  triple  term  will  not  contribute  and  we  find  that  the  corresponding  Sy.^{k)  term  has  the 
limiting  form 


2p(l) 

2 


(58) 


as  k  —*  Q.  Hence 


fi{k  — ♦  0)  = 


_2 _ l_ 

y  +  1  m(i)’ 


(59) 


4.5  The  momentum  equation  eddy  viscosity  and  diffusivity 


The  o.d.e’s,  Eqs.  (45)  and  (48),  for  the  momentum  equation  eddy  viscosity  and  diffusivity  are 
readily  solved  and  shown  in  Fig.  2.  We  observe  that  both  the  eddy  viscosity  and  diffusivity 
have  a  similar  plateau  structure  as  A  — >  0.  Notice  that  the  eddy  diffusivity  plateau  is  not 
obtained  in  the  original  numerical  calculation  of  Rose®.  As  A  — >  Ac,  eddy  viscosity  displays  a 
weak  cusp  like  behavior  while  the  eddy  diffusivity  decreeises  monotonically  as  A  — +  k^.  These 
curves,  as  A  — >  Ac,  are  similar  to  that  of  Zhou  et  al."  and  Rose®. 


5  The  turbulent  transport  coefficient  in  the  second 
moments 


The  concept  of  the  spectral  eddy  viscosity  and  diffusivity  are  introduced  in  the  second 
moments®’'*.  Thus,  the  momentum  equation  spectral  eddy  viscosity  and  diffusivity  are  only 


a  partial  contribution  to  the  total  transport  coeflicienls.  Indeed  from  our  numerira!  mea¬ 
surement  in  Fig.  I,  we  expect  that  the  triple  nonlinear  terms  will  contribute  to  the  energy 
transfer  when  k  is  near  kc- 

We  consider  the  contribution  of  the  triple  nonlinear  term  in  the  renormalized  eddy  vis¬ 
cosity  to  the  eddy  viscosity  first  (Zhou  and  Vahala^^).  The  second  moment  for  the  velocity 
field  is  defined  as 


Uc,0{k,t)  =<  Ua{kj)u0{-kj)  >  . 
The  time  evolution  of  Ua0{k,t)  is 


+  2  <  fjk,t)up{-k,t)  >  +r°(k.<)  +  7-L('<-o  («i) 

In  this  equation,  T^(k,  t)  is  the  standard  energy  transfer  from  the  quadratic  nonlinearity. 
In  contrast,  Tj^(k,  t)  =  —2v'r{k)k'^E{k)  is  the  energy  transfer  arising  from  the  RNG  induced 
triple  nonlinearity.  It  is  readily  shown  that*® 


VT{k)  = 


1  /-fc-tfcc 


/ 

Jkc 


djdz 


Likj^q)\k-i\-y-^jy^^^^ 


(62) 


In  Fig.  3,  we  see  that  i/T(fc)  is  the  major  contributor  to  the  cusp-like  behavior  of  the 
spectral  eddy  viscosity  s,s  k  —*  0.  For  A:  <C  1,  it  has  a  backscatter  of  energy  from  the  subgrid 
and  resolvable  scales. 

We  now  define  Q{k)  as  the  scalar  variance  Q{k)  =<  T{  —  k,t)T{k,t)  >.  The  dynamic 
equation  for  0(fc)  can  be  constructed  by  multiplying  Eq.  (25)  with  T(  — fc,f),  followed  by  an 
averaging  operation.  Again,  a  quasi-normal  approximation  is  applied  to  reduce  the  fourth 
moment  to  the  product  of  the  second  moments.  Notice  that  the  last  term  on  the  RHS  of  Eq. 
(25)  will  not  contribute  to  the  spectral  equation  since  <  u^(j',  t)u.Y(k  —  j  —  j',  < )  >~  ^(k  —  j), 
a  condition  that  can  not  be  satisfied  since  k  —  j  is  in  the  subgrid  scale. 

The  dynamical  equation  for  the  scalar  variance  is 


■^-f-/i(A:)A:^]0(A:,O  =  S^  +  S"'  (63) 

where  is  the  usual  transfer  function  for  the  passive  scalar.  S^{A:)  is  the  additional 
contribution  from  the  triple  nonlinear  term  induced  by  the  recursive  RNG  procedure 


S^(A:) 


2k. 


fi{kc)k 


m-f  1/2 

c 


/  -  jj'-oO  -  j')  X  r(j')r(-k)  > 


2//7’(fc)A:^0(fc,  t)  (64) 


where 


Hrik)  = 


1  f  ,3  .Dcojk  -  j)Q{\k-  j\) 


(65) 


and  the  incompressible  condition  has  been  used.  It  is  seen  in  Fig.  4  that  iirik)  is  small  when 
k  is  small.  However,  as  k  —*  kc,  firik)  increases  rapidly. 

The  solution  of  urik)  is  very  similar  to  that  of  (irik)  as  k  —>■  k^.  They  are  the  major 
contribution  to  the  strong  cusp  in  the  eddy  viscosity  found  from  the  Test  Field  modeP  and 
EDQNM"^.  Furthermore,  i^rik)  also  contains  the  backscatter  of  the  energy  from  the  subgrid 
to  the  resolvable  scale.  This  is  a  major  difference  between  firik)  and  t'rik)- 

Rose*^  discussed  the  role  of  the  triple  nonlinear  terms  in  physical  space.  He  pointed  out 
that  it  represents  the  possibility  of  an  exchange  of  scalar  eddies  between  the  resolvable  and 
subgrid  scales.  This  effect  is  an  inherent  property  of  metisurernents  made  on  the  passive 
scalar  system  with  instruments  which  have  a  spatial  resolution  limited  to  an  eddy  w'idth  size 
greater  than  l/k^. 


6  Spectral  eddy  viscosity,  difFusivity  and  Prandtl  num¬ 
ber 

The  spectral  eddy  viscosity  is  simply  the  sum  of  the  contributions  from  the  momentum 
equation  and  the  triple  nonlinear  term.  The  result  is  presented  in  Fig.  5.  It  appears  that 
our  calulation  is  in  qualitative  agreement  with  the  closure  theory^'"*  and  direct  numerical 
measurements'**-**.  In  particular,  it  predicts  the  correct  asymptotic  behaviors  of  the  eddy 
viscosity  as  k  0  and  k  — »  k^  (Kraichnan*). 

Our  spectral  eddy  diffusivity  shows  a  plateau  at  k  — >  0,  in  good  agreement  with  the 
EDQNM  calculation  of  Chollet**.  However,  the  EDQNM  calculation  is  not  unique  and  de¬ 
pends  on  the  choice  of  the  parameters.  Our  diffusivity  is  in  good  agreement  with  EDQNM 
when  parameters  are  chosen  according  to  the  direct  interaction  approximation  (DIA)**. 

The  spectral  Prandtl  number  can  be  easily  determined  from  our  calculated  eddy  viscosity 
and  diffusivity  (Fig.  5).  It  is  a  function  of  k  and  has  values  ranging  from  0.72  ~  0.92.  Note 
that  our  turbulent  Prandtl  number  at  A;  — >  0  limit  is  very  close  to  that  reported  by  Yakhol 
and  Orszag*  (0.7179).  The  values  of  turbulent  Prandtl  numbers  found  experimentally 
in  the  boundary  layer  are  in  the  range  of  0.6  ~  0.8.  Hinze^*^  and  Tennekes  and  Lumley^' 
pointed  out  that  the  transfer  of  passive  scalar  and  velocity  fields  may  be  equally  effective. 
Thus  the  turbulent  Prandtl  number  is  about  1.  However,  Lesieur  and  Rogallo'*  found  that 
their  spectral  Prandtl  number  only  rose  from  0.2  at  small  k  to  0.8  near  the  cutoff.  Lesieur''* 
recently  found  that  the  turbulent  Prandtl  number  may  be  much  closer  to  1  than  that  of 
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Lesieur  and  Rogallo’*.  The  EDQNM  spectral  Prandtl  number  depends  on  the  choice  of  two 
EDQNM  closure  scalar  adjustable  constants'*’^'’.  In  the  first  case,  it  remains  approximately 
equal  to  0.6,  even  in  the  vicinity  of  kc-  In  the  second  case,  it  has  a  plateau  value  of  0.33, 
and  a  cusp  close  to  Atc,  where  it  rose  to  0.6. 

7  Discussion  and  conclusions 

In  this  paper  we  have  applied  recursive  RNG  to  the  problem  of  the  advection  of  a  passive 
scalar  by  incompressible  turbulence.  We  have  clarified  the  role  of  the  higher-order  RNG 
induced  nonlinearities  and  shown  that;  (a)  The  renormalized  evolution  equations  are  still 
Galilean  invariant  (i.e.,  these  higher-order  noniinearities  do  not  destroy  the  Galilean  invari¬ 
ance  of  the  original  equations).  This  is  an  important  property  that  needs  to  be  preserved  in 
subgrid  modeling,  especially  as  one  proceeds  to  more  complicated  flows  and  boundaries,  (b) 
These  higher-order  nonlinearities  do  not  contribute  to  the  transport  coefficients  as  A:  — >  0. 

Now  the  topical  byproduct  of  the  recursive  RNG  methods  is  a  complicated  integro- 
difference  recursion  relation  to  be  solved  for  the  eddy  transport  coefficients®’*®"".  This 
recursion  relation  is  a  function  of  the  subgrid  shell  thickness  parameter  /.  Here,  we  have 
shown  how  to  pass  to  the  differential  subgrid  shell  thicknef  limit  f  I .  In  this  limit,  we 
recover  an  ordinary  differential  equation  for  the  eddy  coefficients  -  an  equation  that  is  very 
easily  solved. 

The  o.d.e  that  is  derived  in  recursive-RNG  is  fundamentally  different  from  that  derived 
by  e-RNG  techniques.  For  convenience,  we  have  summarized  the  Yahkot-Orszag  derivation 
of  the  €-RNG  o.d.e  in  the  Appendix.  In  e-RNG,  one  is  forced  into  taking  the  ^  >  0  limit*"®, 
and  the  independent  variable  of  the  resulting  o.d.e  is  actually  the  cut-off  wavenumber  kc-  In 
recursive  RNG,  the  independent  variable  is  the  resolvable  scale  wavenumber  k,  0  <  k  <  kci 
with  a  renormalization  transformation  that  permits  kc  to  be  fixed.  There  is  no  renormal¬ 
ization  transformation  made  in  the  Yakhot-Orszag  c-RNG  formulation.  In  the  limit  i  — >  0, 
the  eddy  transport  coefficients  from  both  theories  are  in  very  close  agreement.  This  is  to 
be  expected  since  the  higher-order  recursive  RNG-induced  nonlinearities  — >  0  as  A:  0. 

The  slight  diflference  in  the  eddy  coefficients  (in  the  k  —>  0)  between  the  two  theories  can  be 
attributed  to  the  treatment  of  kc  :  i.e.,  whether  one  performs  RNG  rescaling  transformations 
(recursive  RNG)  or  not.  The  important  effect  of  the  triple  nonlinearities  introduced  by  the 
recursive  RNG  procedure  are  clearly  seen  in  the  second  moment  equations  -  especially  for  a 
resolvable  wavenumber  k  — »  kc- 

The  spectral  eddy  viscosity,  diffusivity  and  Prandtl  number  are  determined  and  we  find 
good  agreement  with  both  closure  theory®"**  and  direct  numerical  simulations*'*"*®. 
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Appendix:  Brief  summary  of  the  Yakhot-Orszag’s  derivation 
of  o.d.e.  for  the  eddy  viscosity 


In  the  e-  RNG,  a  small  parameter  e  is  introduced  into  the  exponent  of  the  forcing  corre¬ 
lation  function,  Eq.  (4), 

t  =  {Al) 

where  d  =  3  is  the  dimensionality®. 

The  distance  interaction®  is  introduced  at  the  outset®,  k  — >  0.  Thus,  the  triple  nonlinear¬ 
ities  do  not  contribute  to  the  recursion  relation,  Eq.  (14).  Furthermore,  the  eddy  viscosity 
is  a  function  of  the  cutoff  wavenumber  kc  and  is  obviously  independent  of  the  resolvable 
scale  wavenumber  0  <  k  <  k^.  To  emphasize  this,  we  shall  follow  Smith  and  Reynolds^^  and 
denote  their  eddy  viscosity  by  i/{k  — »  0,  kc). 

By  variation  of  the  cutoff  wavenumber  k^  Yakhot  and  Orszag®  derived  an  o.d.e  for 

if{k  — >  0,  kc) 


du{k  0,  kc) 
dkc 


=  A‘3v{k  — >  0,  kc)\^{k  — >  0,  kc) 


{A2) 


where 


^3  = 


As  = 


6  —  e 


X{k  —*  0,  kc)  = 


Do 


(^3) 


(27r)3’  30  [u{k-^0,kc)f/^kl^'' 

where  S3  =  47r^  is  the  area  of  a  unit  sphere  in  three-dimensions.  These  equation  should  be 
compared  with  that  Eqs.  (43)-(45)  where  the  resolvable  scale  wavenumber  k  is  the  variable 
in  the  recursion  relation  derived  from  the  recursive  RNG,  and  not  the  cutoff  kc. 

For  a  given  boundary  condition  u{k  —y  0,  fco),  an  analytical  solution  can  be  obtained  from 
Eqs.  (A1-A3). 


u{k  —y  0,  kc)^  —  i>{k  — >  0,  koY 


3A3D0S3  ‘  -  fcp  ‘ 

(27r)^  e 


(A4) 


where  ko  is  typically  in  the  order  of  Kolmogorov  dissipation  wavenumber.  Eq.  (A4)  reduces 
to  the  familiar  inertial  range  form  by  taking  i/{k  —*  0,  fco)  for  kc  k^.  In  that  limit, 

Eq.  (A4)  takes  the  form 


u{k  -y  0,kc) 


T2S3D0 1 1/3  r  3  A3 1 1/3 


I  (27r)3  J 


2c 


(A5) 


where  u{k  —y  0,ko)  is  negler  ed  since  i/(A;  —y  0,kc)  »  v{k  —y  0,ko)  in  the  fully-developed 
turbulence. 


In  order  to  compare  with  the  inertial  range  theory  of  Kraichnan^'^  .  Yakliut  and  Orszag 
set 

kc  =  k-  (,d(i) 

for  any  k  6  [0,^:,-].  From  the  closure  theory  analysis  and  direct  numerical  simulation 
measurements'"^"'*^,  eddy  viscosity  has  rather  distinct  characteristics  in  the  limits  k  0  and 
k  — >  kc-  As  a  result,  Eq.  (A6)  is  a  rather  oversimplified  approximation^^. 
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Fig.  1  Forced  eddy  viscosity  profiles  determined  from  LES  databases  for  the  fluid  velocity 
at  one  time  instant.  where  i/^^{k)  arises  from  measured  LE.S 

nonlocal  subgrid  energy  transfer,  and  arises  from  measured  LES  local  subgrid  energy 

transfer.  It  is  important  to  note  that  i/^“^(k)  — ►  0  as  A:  — +  0  and  that  i/^^{k)  arises  from  the 
u^-u^  interaction.  Note  also  the  cusp  bahvior  in  i/^'^{k)  as  k  — +  Ay. 


v,{k) 


Fig.  3  The  drain  eddy  viscosity  arising  from  the  triple  nonlinearities  in  the 

differential  subgrid  shell  limit  in  recursive  RNG.  r  =  kd Kp  is  a  parameter  in  the  production- 
type  energy  spectrum,  so  that  E{k)  — +  fc”*  as  — >  0.  Kp  is  a  parameter  that  controls 
the  location  of  the  peak  in  E{k).  As  r  increases,  this  peak  in  E{k)  moves  to  smaller  k. 
Backscatter  of  energy  from  the  subgrid  scales  to  the  large  spatial  scales  is  seen  for  k/k^  <  0.4, 
the  region  in  which  urik)  <  0.  For  r  >  1,  there  is  a  sharp  cusp  as  k  — >  k^- 
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r  is  as  in  Fig.  3.  Notice  that  there  is  now  no  backscatter  of  scalar  variance,  since  fiT(k] 
non-negative  for  all  k.  There  is  a  strong  cusp  as  k  k^. 
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